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10Hz  a.c.  current  densities  up  to  1  ma/cm  ,  .chondrocytes  in  hydrated  collagen  gels 
respond  to  current  densities  as  low  as  50  ya/cm^  Hz.  This  suggests  that  either 
the  enzymatic  digestion  process  alters  the  behavior  of  cells  or  the  interaction  of 
the  cells  with  the  extracellular  matrix  modulates  the  response. 

As  we  have  demonstrated  a  decrease  in  protein  synthesis,  one  is  left  to  question  why 
the  cell  reduces  its  protein  synthesis,  and  if  another  cell  function  has  been  enhanced. 
Studies  on  the  synthesis  of  glycosaminoglycans  by  chondrocytes  in  hydrated  gels  suggest 
that  GAG  production  is  also  depressed.  We  have  recently  begun  to  utilize  DNA  staining 
techniques  with  flow  cytometric  measurements  to  investigate  mitotic  activity  in  the 
presence  of  very  low  frequency  electric  fields. 

Other  research  being  examined  is  the  altering  of  protein  secretion  rate  of  cells  by 
modulations  of  exocytosis.  Under  the  guidance  of  Dr.  David  Albertini  (consultant), 
we  have  established  an  experimental  method  and  performed  initial  measurements  of 
the  rate  of  transport  of  membrane  vesicles  within  bovine  fascia  and  human  dermal 
fibroblasts.  The  vesicles  are  labeled  with  Lucifer  Yellow  CH  (1  mg/ml)  which  is 
incorporated  into  the  vesicles  during  serum  stimulated  liquid  phase  endocytosis. 

Thus  far,  only  a  few  experiments  have  been  done  and  no  difference  in  the  rates  of 
transport  of  vesicles  has  been  detected. 

Efforts  are  well  underway  to  establish  a  real-time  fluorescent  assay  of  intracellular 
NAD(P)H  levels.  In  theory,  changes  in  cytoplasmic  NAD(P)  H  levels  should  be  correlated 
to  changes  in  the  rate  of  protein  biosynthesis.  Therefore,  as  proposed  this  technique 
should  enable  us  to  sense  changes  in  cellular  biosynthesis  in  real  time.  With  the 
radiolabeled  precursor  incorporation  studies  as  the  benchmark,  we  plan  to  extend  the 
frequency-amplitude  threshold  response  determinations  to  finely  delineate  the  cutoff 
frequencies  using  this  new  method.  This  should  provide  invaluable  information 
regarding  the  kinetics  of  the  electrochemical  transductive  coupling  at  the  plasma 
membrane . 

Thus  far,  the  installation  of  the  vibration  isolation  table,  pulsed  nitrogen  laser, 
inverted  microscope,  emission  spectrometer,  dual  channel  boxcar  integrator  and  other 
equipment  needed  to  measure  the  emitted  fluorescent  light  from  a  single  cell  excited 
at  340  nm  (near  UV)  wavelength  has  been  completed.  Calibration  using  buffered 
solutions  of  NADH  (Sigma  Chem)  has  been  performed. 

Finally,  efforts  have  been  undertaken  to  determine  the  relationship  between  an 
applied  source  and  the  electric  field  generated  between  two  electrodes.  The 
rectification  and  harmonic  distortion  with  respect  to  a  pure  sinusoid  that  occur 
when  voltage  and  current  sources  are  applied  across  the  electrodes  have  and  are 
being  quantified  and  explained. 


I.  SYNOPSIS 


-Project  objectives  for  the  first  year  include: 

1.  Define  the  relationship  between  physiologic  electrical 
currents  and  the  rate  of  biosynthesis  of  extracellular  matrix 
proteins  by  bovine  fibroblasts  and  chondrocytes. 

2.  Determine  the  extent  of  cellular  proliferative  response  to 
the  same  physiologically  relevant  electrical  stimulus. 

3.  Set  up  and  calibrate  a  real-time  fluorometric  assay  of 
intracellular  NAD(P)H  using  microfluorometric  techniques^ 


I. a  Relevance 

The  long  term  goals  of  this  project  are:  (1)  to  determine  which 
cellular  biochemical  processes  are  modulated  by  physiologic  electric 
fields,  (2)  to  determine  the  physiochemical  transduction  mechanisms  which 
couple  extracellular  fields  to  intracellular  events,  and  (3)  to  explore  the 
usefulness  of  electric  stimuli  as  a  control  of  synthetic  tissue  growth  and 
development . 

Our  primary  working  hypothesis  is  that  many  mammalian  non-excitable 
cell  types  are  sensitive  to  naturally  occurring  interstitial  electrical 
currents  which  arise  from  in-vivo  sources  and  that  cells  use  these  signals 
to  establish  and  maintain  the  structural  integrity  of  their  environment. 
Evidence  supporting  this  concept  is  substantial.  Furthermore,  we  find  it 
plausible  that  electric  stimuli  may  be  useful  to  regulate  generation  of 
living  tissue  substitutes,  modulation  of  normal  wound  healing  and  perhaps 
in  the  management  of  hypertropic  disorders  of  wourtd  healing.  We  believe 
these  goals  to  be  consistent  with  the  mission^of  the  Office  of  Naval 
Research.  -  •" 


I.b  Summary 


>With  an  in-vitro  system  which  minimized  the  sample-to-sample 
variation  by  using  a  isogenetic  cell  population  at  a  uniform  cell  density, 
we  observed  that: 


Bovine  (and  Human)  fibroblasts  and  bovine  chondrocytes  respond  to  low 
frequency  sinusoidal  electric  current  by  altering  the  rate  of  biosynthesis 
of  extracellular  matrix  structural  macromolecules.  This  effect  was 
dependent  on  stimulus  frequency  and  amplitude. 


52^&'The  sensitivity  of  the  fibroblasts  to  the  fields  was  found  to  depend  on 
the  orientation  of  the  cells  in  the  direction  of  the  field.  Cells  oriented 
such  that  their  major  axis  was  aligned  in  the  direction  of  the  electric 
field  responded  to  weaker  electric  currents  than  those  aligned 
perpendicular  to  the  field.  This  evidence  is  consistent  with  the  hypothesis 
that  the  mechanism  of  electrochemical  transduction  involves  a  change  in  the 
transmembrane  potential. 


r~~ 

Jil  We  have  also  observed  that  while  chondrocytes  in  organ  culture  do  not 
change  their  biosynthetic  behavior  in  response  to  10  Hz  a.c.  current 
densities  up  to  1  ma/dm2,  chondrocytes  in  hydrated  collagen  gels  respond  to 
current  densities  as  low  as  50  pa/cm2  at  20  Hz*  This  suggests  that  either 
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1  the  enzymatic  digestion  process  alters  the  behavior  of  cells  or  the 
interaction  of  the  cells  with  the  extracellular  matrix  modulates  the 
response. 

As  we  have  demonstrated  a  decrease  in  protein  synthesis,  one  is  left 
to  question  why  the  cell  reduces  its  protein  synthesis,  and  if  another  cell 
function  has  been  enhanced.  Studies  on  the  synthesis  of  glycosaminoglycans 
by  chondrocytes  in  hydrated  gels  suggest  that  GAG  production  is  also 

depressed.  We  have  recently  begun  to  utilize  DNA  staining  techniques  with 
flow  cytometric  measurements  to  investigate  mitotic  activity  in  the 

presence  of  very  low  frequency  electric  fields. 

Other  research  being  examined  is  the  altering  of  protein  secretion 
rate  of  cells  by  modulations  of  exocytosis.  Under  the  guidance  of  Dr.  David 
Albertini  (consultant),  we  have  established  an  experimental  method  and 
performed  initial  measurements  of  the  rate  of  transport  of  membrane 

vesicles  within  bovine  fascia  and  human  dermal  fibroblasts.  The  vesicles 

are  labeled  with  Lucifer  Yellow  CH  (1  mg/ml)  which  is  incorporated  into  the 
vesicles  during  serum  stimulated  liquid  phase  endocytosis.  Thus  far,  only  a 
few  experiments  have  been  done  and  no  difference  in  the  rates  of  transport 
of  vesicles  has  been  detected. 


Efforts  are  well  underway  to  establish  a  real-time  fluorescent  assay 
of  intracellular  NAD( P)H  levels.  In  theory,  changes  in  cytoplasmic  NAD(P)H 
levels  should  be  correlated  to  changes  in  the  rate  of  protein  biosynthesis. 
Therefore,  as  proposed  this  technique  should  enable  us  to  sense  changes  in 
cellular  biosynthesis  in  real  time.  With  the  radiolabeled  precursor 
incorporation  studies  as  the  benchmark,  we  plan  to  extend  the 
frequency-amplitude  threshold  response  determinations  to  finely  delineate 
the  cutoff  frequencies  using  this  new  method.  This  should  provide 
invaluable  information  regarding  the  kinetics  of  the  electrochemical 
transduct ive  coupling  at  the  plasma  membrane. 

Thus  far,  the  installation  of  the  vibration  isolation  table,  pulsed 
nitrogen  laser,  inverted  microscope,  emission  spectrometer,  dual  channel 
boxcar  integrator  and  other  equipment  needed  to  measure  the  emitted 
fluorescent  light  from  a  single  cell  excited  at  3^0  nm  (near  UV)  wavelength 
has  been  completed.  Calibration  using  buffered  solutions  of  NADH  (Sigma 
Chem.)  has  been  performed. 

Finally,  efforts  have  been  undertaken  to  determine  the  relationship 
between  an  applied  source  and  the  electric  field  generated  between  two 
electrodes.  The  rectification  and  harmonic  distortion  with  respect  to  a 
pure  sinusoid  that  occur  when  voltage  and  current  sources  are  applied 
across  the  electrodes  have  and  are  being  quantified  and  explained. 


II.  BIOSYNTHETIC  RESPONSES 


Introduction 


The  role  of  physiologic  electric  fields  in  the  regulation  of  tissue 
growth,  remodeling,  and  repair  present  questions  of  burgeoning  interest  to 
biologists  and  clinicians.  We  studied  the  sensitivity  of  mammalian 
fibroblasts  to  electric  currents  using  neonatal  bovine  fibroblasts  in 
collagen  gels.  The  frequency  and  intensity  dependence  of  electric  field 
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modulation  of  proline  incorporation  into  extracellular  and  intracellular 
protein  was  measured  over  the  frequency  range .  of  0.1  Hz  to  1  kHz.  A 
frequency  and  amplitude  dependent  reduction  in  the  rate  of  incorporation 
was  observed.  Using  tissues  containing  cells  aligned  either  parallel  or 
perpendicular  to  the  electric  field,  it  was  observed  that  this  response  was 
dependent  on  the  orientation  of  the  cells  relative  to  the  direction  of  the 
applied  electric  field. 

Determination  of  cell  sensitivity  to  electric  fields  which  span  the 
physiologic  frequency  range  would  allow  an  assessment  of  the  extent  to 
which  physiologic  fields  could  regulate  cellular  processes.  Comparison  of 
frequency  dependence  with  the  kinetics  of  known  cellular  processes  may 
provide  insight  to  the  putative  mechanisms  of  energy  transduction. 

II. b  Methods 

Connective  tissue  remodeling  is  known  to  be  regulated  by  physical 
stresses  [1],  Since  fibroblasts  are  primarily  responsible  for  this 
remodeling  activity,  these  cells  were  selected  for  use  in  this  study. 
Neonatal  bovine  fibroblasts  were  obtained  by  disaggregating  superficial 
fascial  tissue  from  the  thigh  of  2  week  old  calves  by  serial  trypsin  and 
collagenase  digestions.  The  cells  obtained  were  plated  in  flasks, 
maintained  in  Dulbecco's  Modified  Eagle's  medium  (DME)  containing  10$  calf 
serum  and  passed  2  to  5  times  prior  to  incorporation  into  gels. 

The  behavior  of  many,  if  not  all  cell  types,  is  strongly  influenced 
by  the  cells'  environment.  So,  to  maximize  control  over  extracellular 
matrix  composition  and  cell  density,  we  chose  to  fabricate  tissues  of 
constant  composition  by  incorporating  fibroblasts  in  collagen  matrices 
( FPCM) .  To  fabricate  the  fibroblast  populated  collagen  matrices  (FPCM), 
bovine  fibroblasts  were  incorporated  into  2  mg/ml  type  I  collagen  gels 
using  the  technique  reported  by  Bell  et  al.  [2].  Native  collagen  was 
obtained  through  0.5M  acetic  acid  extraction  of  the  tail  tendons  of  young 
(less  than  2  month  old)  Sprague-Dawley  rats.  The  soluble  collagen  solution 

was  stored  in  sterile  0.1M  acetic  acid  at  4  °C  until  needed.  At  the  time  of 
gel  casting  neutralized  collagen  solution  and  cells  suspended  in  media  were 

mixed  and  warmed  to  37  °C.  Gelation  occurred  within  10  minutes  and  within 
several  hours  the  gel  was  visibly  contracting.  Three  days  after  gel  casting 
both  matrix  contraction  and  cell  population  doubling  were  found  to  plateau 
and  the  FPCM  were  strong  enough  to  be  handled.  Hematoxylin  and  eosin  (H&E) 
staining  of  10  micron  thick  histologic  sections  [  McLeod,  submitted  ] 
revealed  a  randomly  oriented  FPCM  with  cell  and  matrix  densities  similar  to 
that  of  loose  connective  tissue.  Cells  in  these  FPCM  stopped  proliferating 
prior  to  contact  inhibition  and  cell  density  approximated  loose  connective 
tissue  [2],  The  electrical  conductivity  of  the  FPCM  was  found  to  be  less 
than  10$  different  from  the  surrounding  media.  On  the  third  day  after 
casting,  1  cm2  samples  were  cut  from  the  FPCM  and  randomly  assigned  to 
control  or  exposure  groups. 

Control  and  exposure  samples  were  mounted  in  teflon  holders  which  had 
circular  apertures,  and  then,  mounted  in  a  teflon  chamber  (Fig.  1) 
containing  fresh  serum-free  media  supplemented  with  10  yCi/ml  3H-proline 
and  1  millimolar  proline.  The  chambers  were  incubated  in  95$  air-5$  CO2 

mixture  at  37  °C  and  99$  humidity.  A  programmable  current  source  provided 
stable  low  frequency  currents  with  d.c.  flow  limited  to  less  than  0.1$  of 
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the  a.c .  amplitude.  Current  was  passed  through  the  exposure  samples  for  12 
hours  via  platinum  electrodes  which  were  separated  from  the  bath  by  media 
bridges  and  convection  barriers  (membranes).  Current  density  was  determined 
by  dividing  the  total  current  by  the  aperture  area  of  the  holders,  and 
field  intensity  in  the  FPCM  samples  was  estimated  using  Ohm's  Law  and  the 
value  of  the  media  resistivity. 


- 


Figure  1  :  Experimental  apparatus  used  to  expose  random  FPCM  samples  to 
electrical  fields.  Platinum  electrodes  are  separated  from  the  site  of  the 
tissue  sample  via  media  bridges.  A  programmable  current  source  with  series 
d.c .  blocking  capacitors  maintains  a  constant  current  within  the  chamber 
media,  independent  of  variations  in  electrode-electrolyte  interfacial 
impedance.  The  lumen  of  the  sample  holder  establishes  a  defined  current 
density,  and  therefore  field  intensity,  through  the  FPCM 


Accumulation  of  newly  synthesized  Drotein  was  estimated  by  measuring 
the  incorporation  of  proline  into  macror olecules .  Immediately  after  a  12 

hour  exposure  period,  samples  were  washed  at  4  °C,  then  digested  with 
bacterial  collagenase,  after  which  the  cells  were  separated  from  the  matrix 
components  by  centrifugation.  After  counting  the  cells,  cell  were  disrupted 
[31  and  intracellular  protein  precipitated  with  cold  trichloroacetic  acid. 
Scintillation  counts  of  both  the  matrix  and  cell  fractions  were  obtained 
and  normalized  to  cell  number.  Control  studies  demonstrated  an 
incorporation  rate  of  450  picomoles  or  proline  per  hour  per  million  cells 
in  serum-free  medium,  addition  of  10$  calf  serum  increased  this  synthesis 
rate  by  a  factor  of  1. 7-2.0.  Collagen,  assayed  by  differential  salt 
precipitation  [4],  accounted  for  approximately  4$  of  recovered  protein 
synthesized  [5]  which  is  consistent  with  a  previous  report  [6]. 

Experiments  were  performed  over  a  range  of  current  densities  (  0.1 
ya/cm2  to  1  ma/cm2,  r.m.s.)  and  frequencies  (  0.1  to  1000  Hz  ).  The 
observed  response  to  electrical  current  was  an  alteration  in  the  rate  of 
incorporation  of.  3H-proline  into  secreted  protein.  Figure  2a  exhibits  the 
response  as  a  function  of  current  density  at  1  Hz  and  shows  that  currents 
as  low  as  1  yA/cm2  cause  a  30$  reduction  in  the  incorporation  of  newly 
synthesized  protein  into  the  extracellular  matrix.  Above  this  current 
density,  increasing  the  current  over  two  orders  of  magnitude  does  not 
significantly  increase  the  response.  However  halving  the  current  density 
leads  to  no  response.  The  transition  is  described  as  a  threshold  because 
the  response  is  strongly  field  intensity  dependent  for  only  a  very  small 
fraction  of  the  total  range  of  intensities  which  we  found  capable  of 
evoking  the  response.  The  threshold  current  density  was  found  to  be 
frequency  specific. 


The  frequency  dependence  of  the  response  plotted  in  Figure  2(b) 
suggests  a  frequency  range  which  may  be  optimally  effective  for  alteration 
of  extracellular  matrix  protein  synthesis.  This  observation  is  summarized 
in  Figure  2c  where  the  effect  of  both  frequency  and  current  density  are 
plotted  over  the  range  of  1-1000  Hz  and  0.1  to  100  yA/cm2  respectively. 
Peak  sensitivity  was  recorded  at  10  Hz  where  a  current  density  of  only  0.5 
yA/cm2  produced  a  significant  reduction  in  protein  incorporation  in  the 
matrix  component.  Since  the  media  resistivity  was  65  ft-cm  [7],  this 
corresponds  to  a  peak  field  intensity  in  the  medium  of  45  yV/cm.  The 
incorporation  of  radiolabel  into  intracellular  protein  was  found  to  reflect 
the  pattern  seen  in  the  extracellular  matrix. 


To  test  the  possibility  that  chemical  by-products  of  electrolysis 
were  responsible  for  the  response,  experiments  were  performed  in  which  100 
ya/cm2  was  passed  through  the  media  for  12  hours,  then  discontinued 
immediately  prior  to  transfer  of  the  FPCM  to  the  chambers.  The  FPCM  were 
then  kept  in  the  chambers  for  12  hours  while  bathed  in  the  electrically 
conditioned  media  for  12  hours  then  analyzed  as  usual.  No  significant 
difference  in  rates  of  proline  incorporation  was  observed  between  the 
control  samples  and  samples  which  had  been  in  the  electrically 
preconditioned  media.  This  suggested  that  the  responses  measured  required 
direct  exposure  to  the  electric  field  and  the  electric  field  effects  were 
not  mediated  by  stable  non-volatile  by-products  of  electrolysis.  Similarly, 
to  test  whether  the  cellular  responses  were  a  consequence  of  thermal 
effects,  the  magnitude  of  the  Joule  heating  was  analyzed.  The  temperature 
rise  of  the  media  adjacent  to  the  sample  position  was  measured  as  a 
function  of  both  current  density  and  time.  It  was  found  that  a  current 


density  of  10  ma/cra2  was  required  to  produce  a  steady  state  temperature 

increase  of  0.5°C  with  an  onset  time  constant  approximately  equal  to  20 
minutes  [7].  The  steady  state  temperature  elevation  was  proportional  to  the 
current  density  squared.  For  the  range  of  current  densities  used  in  the 

actual  experiments,  the  expected  temperature  rise  was  much  less  than  0.1°C, 
which  is  less  than  the  temperature  oscillations  in  the  incubators. 
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Figure  2:  Response  of  bovine  fibroblasts  in  FPCM  to  12  hour  exposures  to 


low  frequency  electric  fields,  at  37  C,  in  serumless  Dulbecco's  media. 
3H-proline  incorporation  (cpm)  were  normalized  to  cell  count  for  each 
sample.  The  incorporation  ratio  of  experimental  to  control  samples  was 
determined  and  the  natural  logarithm  (Ln)  is  statistically  analyzed.  The 
mean  and  standard  error  of  the  mean  for  n  samples  (6  <  n  <  20  )  are 
plotted,  "p"  value  indicates  probability  of  log  normalized  data  differing 
from  0.0  as  determined  by  the  non-parametric  Dixon  -  Mood  sign  test,  (a) 
Extracellular  proline  accumulation  vs.  current  density  at  1  Hz.  stimulus 
frequency,  (b)  Normalized  extracellular  proline  accumulation  vs.  frequency 
for  1  uA/cm2.  (c)  Minimum  field  intensity  for  a  detectable  response. 
Summary  of  results  for  all  tested  frequencies  and  current  densities. 
Current  densities  converted  to  peak  field  intensities  using  the  measured 
media  resistivity  of  65  Q-cm.  The  lower  boundary  of  the  gray  region 
represents  the  highest  field  intensity  at  which  no  significant  change  in 
extracellular  protein  accumulation  was  detected,  the  upper  boundary 
represents  the  lowest  intensity  evoking  a  significant  change. 


II. c  Dependence  on  Cell  Orientation 

For  the  frequencies  used  in  these  experiments,  plasma  membrane 
impedance  is  many  orders  of  magnitude  higher  than  that  of  the  cytoplasm. 
When  a  cell  the  size  of  a  bovine  fibroblast  is  placed  in  an  applied  field, 
the  membrane  essentially  prevents  current  flow  through  the  cell.  Current 
exclusion  results  in  an  additional  voltage  drop  across  the  cell  membrane 
which  is  spatially  distributed  over  the  cell  surface  [8].  This  distribution 
is  a  function  of  cell  geometry,  such  that  for  any  point  on  the  membrane  the 
induced  perturbation  of  the  membrane  potential  is  proportional  to  EqL/2, 

where  Eq  is  the  applied  field  magnitude,  and  L  the  maximum  linear  dimension 

of  the  cell  in  the  direction  of  the  field.  For  non-spherical  cells,  a 
maximum  perturbation  of  the  membrane  potential  will  occur  when  the  major 
axis  of  the  cell  is  aligned  in  the  direction  of  the  imposed  electric  field. 
The  importance  of  the  cell  membrane  potential  in  regulating  cellular 
activity  suggests  that  this  alteration  in  transmembrane  potential  may 
mediate  the  cellular  response  to  low  frequency  extracellular  electric 
fields . 

One  test  of  this  hypothesis  makes  use  of  the  nonspherical  shape  of 
the  cells  used  in  this  study.  Fibroblasts,  in  the  collagen  lattices,  took 
on  the  bipolar  morphology  which  has  been  previously  described  [9].  The 
maximum  length  of  these  cells  was  approximately  150  pm,  which  is  7-10  times 
longer  than  their  minor  dimension.  If  the  depression  in  incorporation  is 
mediated  by  a  membrane  potential  change,  then  cells  exposed  to  fields 
parallel  to  their  major  axis  should  exhibit  a  lower  threshold  intensity 
than  cells  with  their  major  axis  perpendicular  to  the  applied  electric 
field.  Of  course  other,  probably  significant,  parameter  changes  occur  with 
changes  in  the  cells'  orientation.  The  plasma  membrane  area  over  which  the 
maximum  transmembrane  potential  alteration  occurs  is  reduced  when  cells  are 
aligned  with  the  field.  Also,  the  interaction  between  the  cell  surface  and 
the  tangent  electric  field  will  depend  on  cell  orientation.  Therefore, 
interpretation  of  orientation  dependence  may  be  complicated.  Nonetheless, 
such  an  effect  would  provide  further  evidence  that  the  response  is 
dependent  only  on  the  electric  field. 

We  investigated  the  role  of  cell  orientation  by  constructing  FPCM 
with  cells  predominantly  oriented  in  one  direction.  To  uniformly  align 
cells  in  the  collager  matrix,  the  FPCMs  were  allowed  to  contract  over  3 
days  around  two  porous  polyethylene  posts  held  at  a  fixed  separation 
distance  of  two  centimeters  (Figure  3)*  The  FPCM  contracted  with  the  cells 
and  collagen  aligned  along  an  axis  defined  by  the  line  passing  through  the 
porous  posts  at  opposite  ends  of  the  petri  dish  (Figure  4).  The  electrical 
conductivity  of  the  media  filled  posts  were  within  10$  of  free  media 
solution. 


Figure  3=  Contraction  nad  formation  of  oriented  FPCM.  Fibroblast  populated 
gels  cast  into  60  mm  petri  dishes  with  two  polyethylene  posts  held  at  a 
fixed  separation,  (a)  Hydrated  gel  after  12  hours  of  contraction;  (b)  2U 
hours;  (c)  U8  hours;  and  (d)  72  hours,  the  condition  at  which  the  oriented 
FPCMs  were  utilized. 


Figure  Bovine  fibroblasts  in  oriented  FPCM  as  seen  by  bright-field 
microscopy  18  hours  post  casting.  Vertical  axis  of  photo  is  longitudinal 
axis  between  fixed  Dolvethvlene  Dosts  of  orientation  aDDaratus. 


Three  days  after  casting,  the  FPCM  were  placed  in  the  exposure 
chamber  shown  in  Figure  5.  In  each  chamber,  half  of  the  FPCM  were  installed 
with  the  major  axis  of  the  cells  parallel  to  the  applied  electric  field, 
while  the  other  half  were  installed  with  cells  oriented  perpendicular  to 
the  direction  of  the  electric  field.  Current  was  passed  through  the 
experimental  samples  for  12  hours.  Following  the  exposure  period,  the 
center  section  between  the  posts  were  removed  and  analyzed  for  proline 
incorporation  using  the  same  protocol  as  used  for  the  FPCM  with  random  cell 
orientation. 

Since  the  cells  were  most  sensitive  to  10  Hz  fields,  this  frequency 
was  chosen  to  examine  the  effect  of  orientation  on  the  field  intensity 
threshold.  For  the  randomly  oriented  FPCM,  a  current  density  of  0.3  uA/cm2 
produced  no  significant  effect  on  the  rate  of  proline  incorporation. 
Whereas,  when  cells  aligned  with  the  electric  field  were  exposed  to  the 
same  current  density  a  significant  reduction  in  proline  incorporation  was 
detected.  By  contrast,  the  cells  oriented  perpendicular  to  the  field  did 
not  respond.  As  illustrated  in  Figure  6,  cell  alignment  with  respect  to  the 
electric  field  modulates  the  intensity  threshold.  Cells  parallel  to  the 
field  respond  at  0.3  uA/cm2,  while  cells  perpendicular  to  the  field  have 
still  not  shown  a  significant  depression  in  protein  secretion  at  0.6 
uA/cm2.  At  10  Hz,  the  cells  with  their  major  axis  aligned  with  the  field 
detect  a  field  intensity  of  45  uV/cm.  This  corresponds  to  a  maximum 
membrane  potential  perturbation  of  less  than  0.5  pV. 

II. d  Discussion  of  Biosynthetic  studies 

This  study  has  demonstrated  that  protein  production  in  FPCM  is  more 
sensitive  to  electric  fields  over  the  physiologic  frequency  range  than 
previously  shown  for  connective  tissue  cells.  If  the  same  level  of 
sensitivity  exists  in-vivo  it  seems  possible  that  mechanically  induced 
fields  in  connective  tissue  are  capable  of  triggering  cell  mediated  changes 
in  tissue  composition.  Similarly,  fibroblasts  in-vi vo  may  be  sensitive  to 
electric  fields  from  other  sources,  depending  on  its  intensity  and 
frequency. 

These  experiments  further  demonstrate  that  the  orientation  of  the 
cells  with  respect  to  the  imposed  electric  field  affects  the  ability  of  the 
cells  to  respond.  This  result  is  consistent  with  the  possibility  of  the 
cells  detecting  a  transmembrane  potential  perturbation,  but  it  is  not 
conclusive.  Since  the  FPCM  oriented  parallel  and  perpendicular  to  the  field 
were  in  the  same  bathing  media,  these  studies  strongly  suggest  that  the 
response  was  due  exclusively  to  the  imposed  electric  fields,  and  not  a 
secondary  factor  such  as  heating,  media  contamination,  mechanical 
vibration,  or  ambient  magnetic  fields  within  the  incubators. 

The  frequency  dependence  of  the  response  suggests  two  modes  of 
electrically  mediated  control  of  tissue  composition.  At  constant  frequency, 
a  biosynthetic  response  could  follow  an  increase  in  local  current  density. 
Alternately,  if  a  change  in  loading  rate  leads  to  a  change  in  the 
frequency,  then  this  could  also  lead  to  a  response.  Through  either  pathway, 
the  fields  might  be  used  as  a  feedback  signal  for  tissue  remodeling  and 
repair  [10]. 
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Figure  5:  Exposure  chamber  for  oriented  FPCM  samples.  Each  chamber  holds 


four  samples  with  cells  oriented  parallel  with  the  applied  field,  and  four 
perpend:  ~  *-o  the  applied  field.  Current  density  is  dictates  by  media 

depth  f."  i  -  stant  current  source. 
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Figure  6:  Normalized  extracellular  proline  incorporation  (  closed  symbols  ) 
and  intracellular  proline  incorporation  (open  symbols  )  vs.  current  density 
for  oriented  FPCMs  at  10  Hz.  Samples  with  cells  oriented  parallel  to  field 
demonstrate  a  threshold  current  density  below  that  seen  in  the  samples  of 
randomly  oriented  cells,  while  samples  with  cells  perpendicular  show  no 
depression  of  proline  incorporation  at  a  current  intensity  above  the 
threshold  level  for  randomly  oriented  samples.  Data  points  are  the  mean  and 
standard  error  of  the  mean. 


The  molecular  mechanisms  through  which  weak  electric  fields  trigger  a 
biosynthetic  response  are  unknown.  Since  a  variety  of  cell  types,  including 
the  special  electroreceptors  of  aquatic  animals  [11],  also  exhibit  extreme 
sensitivity  to  electric  fields,  we  suggest  that  this  capability  is  very 
primative  and  serves  a  fundamental  role  in  the  interaction  of  living 
systems  with  their  environment. 


III.  PINOCYTOTIC  VESICLE  TRANSPORT 

III. a  Introduction 

We  have  recently  demonstrated  that  physiologically  relevant  electric 
fields  (those  of  low  magnitude  and  frequency)  affected  the  ability  of 
bovine  fibroblasts  to  incorporate  3H-proline  into  protein.  Secretion  is  the 
final  step  in  the  production  of  extracellular  macromolecules.  Since 
secretion  of  macromolecules  is  partly  mediated  by  exocytosis,  it  is 
plausible  that  electric  fields  might  affect  the  rate  of  biosynthesis 
through  the  alteration  of  intracellular  transport  of  membrane  vesicles. 

To  explore  this  possibility  we  have  undertaken  a  three  phase  project. 
The  first  part  was  the  development  of  a  real  time  technique  for  measuring 
membrane  transport.  The  second  phase  was  developing  a  theoretical  framework 
for  interpreting  membrane  transport  data.  These  have  been  completed.  The 
final  stage,  which  is  currently  being  pursued,  is  the  use  of  the  real  time 
technique,  and  the  use  of  fixed  cells  to  determine  the  potential 
relationship  between  a  field  and  transport.  In  this  report  we  will 
summarize  the  conclusions  of  the  first  two  stages  of  this  effort. 

Ill.b  Background 

Traditionally,  the  three  membrane  transport  processes,  intracellular, 
exocytosis,  and  endocytosis,  are  monitored  using  fixed  cells.  The  cells  are 
allowed  to  utilize  labeled  vesicles  for  various  amounts  of  time  and  then 
are  fixed.  After  fixation  all  of  the  measurements  are  made. 

Fixation  has  many  nice  features.  Sterile  technique  is  not  required 
after  fixation.  Because  the  cells  do  not  have  to  observed  during  the 
experiment,  a  large  number  of  cells  can  be  tested  at  one  time,  and  because 
it  is  so  widely  used,  it  is  the  most  readily  comparable  form  of  data. 

What  data  from  fixed  cells  lack,  though,  is  the  insight  a  real  time 
approach  gives.  In  real  time  an  individual  cell  is  traced  through  time;  it 
is  monitored  while  it  is  alive.  So  a  visual  picture  of  evolution  over  time 
of  an  effect  on  a  single  cell  is  developed.  To  get  such  a  picture  with 
fixation,  one  must  superimpose  the  reactions  of  many  different  cells.  This 
leaves  the  distinct  possibility  that  the  general  idea  might  get  lost  in  the 
noise  of  individual  cellular  differences.  For  this  reason,  it  was  decided 
to  develop  a  real-time  technique  to  measure  the  effect  of  an  electric  field 
on  membrane  transport. 

III.c  Method 

The  protocol  for  tracking  membrane  transport  in  real  time  was  as 
follows: 


Bovine  or  Human  Fibroblasts  (depending  on  the  choice  of  experiment)  of 
passage  number  five  or  less  were  plated  onto  sterile  microscope  cover  glass 
(25mm  by  25mm)  and  maintained  in  Dulbecco's  Modified  Media  (Gibco)  with  10$ 
calf  serum.  The  cells  were  kept  in  incubators  at  37oC  and  with  5$  CO 

Twenty- four  hours  before  the  experiment  was  to  take  place  the  media  was 
changed  to  D.M.E.M.  without  calf  serum.  The  cells  were  labeled  for  forty 
minutes  with  1  mg/ml  Lucifer  Yellow  CH  (Molecular  Probes).  The  cells  were 
washed  five  times  with  Hanks  w/out  Calcium  and  Magnesium  (Gibco)  and  were 
then  given  a  final  wash  with  Hanks  w/  Calcium  and  Magnesium  (Gibco)  to 
remove  any  phenol  red  that  attached.  After  the  washing,  the  cover  glass 
with  the  fibroblasts  was  placed  over  a  1.5cra  diameter  hole  in  the  custom 
built  chamber  shown  in  Figure  7.  The  glass  slip  was  affixed  to  the  teflon 
chamber  by  a  thin  layer  of  Bacitracin  Ointment  U.S.P.  (fougera).  The 
chamber  was  then  filled  with  Hanks  with  Calcium  and  Magnesium. 

Transport  was  observed  by  putting  the  chamber  on  the  stage  of  a  Nikon 
Diaphot  microscope.  The  dye  was  excited  by  a  100W  Osram  Mercury  Arc  Lamp 
that  was  filtered  by  the  Nikon  "B"  filter  block.  A  Nikon  FT  35mm  camera  was 
attached  to  the  microscope  and  pictures  were  taken  in  15  minute  intervals. 
A  neutral  density  filter  (90$  intensity  attenuation)  was  placed  in  front  of 
the  arc  lamp,  and  that  light  path  was  shut  off  between  pictures  to  minimize 
bleaching. 


Proper  temperature  was  maintained  by  the  Nikon  Incubator  NP-2  attached 
to  the  microscope  stage.  Dehydration  was  avoided  by  using  tissue  culture 
plastic  to  cover  the  chamber.  The  electric  field  was  applied  by  using  an 
I.E.C.  F34  function  generator  with  a  Kepco  amplifier.  A  1 00n  resistor  was 
attached  between  the  generator  and  the  chamber.  The  voltage  across  the 
resistor  was  then  measured  with  a  Tektronix  2215  60MHz  oscilloscope,  and 
current  was  determined.  The  field  intensity  was  then  calculated. 

Ten  minutes  after  labeling  was  finished,  the  electric  field  was  turned 
on,  and  transport  was  monitored  for  two  hours.  Under  the  same  conditions  a 
control  experiment  was  also  run.  Pictures  were  taken  ten  minutes  after 
labeling  for  two  hours  in  fifteen  minute  intervals. 

Ill.d  A  Theoretical  Framework  for  Interpreting  Membrane  Transport  Data 


To  analyze  any  changes  that  are  seen  when  monitoring  membrane 
transport,  one  must  understand  the  potential  mechanisms  by  which  the  field 
could  affect  a  cell.  There  are  three  main  possibilities: 


1)  The  field's  impact  could  be  indirect.  By  changing  the  transmembrane 
potential  and  surface  charge  distribution,  the  transport  of  ions  across  the 
membrane  could  be  affected.  This  could  trigger  a  biochemical  reaction  that 
might  alter  rate  of  membrane  pinching  and  fusion  or  the  rate  of  vesicle 
transport . 


2)  The  field  could  have  a  direct  effect.  Pinocytotic  vesicles  have  a 
fixed  surface  charge  and  therefore  will  electrophoresis  in  an  electric 
field.  This  could  alter  transport,  since  the  distance  between  the  vesicles 
and  the  cell  membrane  influences  the  probability  of  fusion. 


3)  Or  it  could  be  a  combination  of  the  above. 


To  understand  these  possible  mechanisms,  and  which  one  is  the 
dominant  cause,  an  expression  for  the  electrophoretic  velocity  of  a  vesicle 


inside  a  cell  was  needed. 


Ill.e  Derivation  of  an  Electrophoretic  Velocity  Equation 

Steady  state  electrophoretic  movement  is  due  to  the  balance  between 
viscous  and  electrical  forces.  This  results  in  the  velocity  of  a  particle 
being  proportional  to  the  strength  of  the  imposed  electric  field. 

Neglecting  the  redistribution  of  the  fixed  membrane  charges  due  to  the 
electric  field,  there  are  two  characteristics  of  pinocytotic  vesicles  that 
can  be  used  to  arrive  at  an  electrophoretic  velocity  expression.  First, 
since  the  double  layer  for  a  pinocytotic  vesicle  is  estimated  at  Inm,  while 
its  diameter  ranges  from  ^Onm  -  500nm  (depending  on  the  cell  type),  it  can 
be  concluded  that  the  electrical  double  layer  is  very  thin  compared  to  the 
dimensions  of  our  colloid  (pinocytotic  vesicle).  Second,  since  the 
pinocytotic  vesicle's  membrane  has  a  much  lower  conductivity  than  that  of 
the  cytoplasm,  the  colloid  can  be  considered  insulating.  Since  its  surface 
conductance  is  low,  the  distribution  of  the  intracellular  electric  field 
is,  therefore,  practically  unaffected  by  the  vesicle. 

Under  those  conditions,  neglecting  relaxation,  electrophoretic 
velocity  can  be  expressed  as  [131: 


where  n  is  the  viscosity  coefficient  of  the  intracellular  fluid.  This 
constant  is  unknown,  but  it  is  known  that  it  must  be  less  than  or 

-3 

equal  to  that  of  water  which  is  1.002  *10  kg/sec-m 

e  is  the  cytoplasm's  dielectric  constant,  again  approximately  that  of 

water  —  80e  —  7.08  •  10  ^  (F/m) 

o 

5  is  the  zeta  potential,  which  is  defined  as  the  electrical  potential 
in  the  slipping  plane  between  the  fixed  and  flowing  liquid  of  the 
planar  model  used  to  derive  our  expression.  It  is  approximately  26mv. 

Ec  is  the  cytoplasmic  field. 

Therefore , 


V  =  1.46  •  10_9(As2/kg)  •  E 
e  c 

The  problem  then  is  to  determine  Ec> 

Ill.f  Determining  The  Cytoplasmic  Electric  Field  Using  a  DC  Spherical  Model 

A  DC  spherical  model  can  often  provide  insight  into  the  magnitude  of 
the  intracellular  fields  and  membrane  voltage  drops  induced  by  a  low 
frequency  stimulus.  For  the  conditions  of  this  project  such  a  model,  shown 
in  Figure  8,  would  contain  the  following  parameters: 

R  —  the  radius  at  which  the  cytoplasm  and  the  membrane  meet  —  30um 

[U] 


J 


4 


AT 


■S3 


'(J -i 

t'v 

$ 

:51 


V, 

I 


$ 

$ 

f 


V 

5f 


»> 

‘.%*l 

1 


R  —  the  radius  at  which  the  membrane  and  the  external  media  touch  — 
e 

30.01 pm 

5  —  the  membrane  thickness  (R^  -  R  )  —  lOnm 

e  c 

o  ,  o  ,  a  —  the  conductivities  of  the  cytoplasm,  membrane,  and 
c  rn  e 

“8  -6 

external  media  —  1  ,  10  to  10  ,  and  1  mho/m,  respectively 

e  ,  e  ,  e  —  the  dielectric  constants  of  the  cytoplasm,  membrane,  and 
c  m  e 

external  media  —  7.08  •  10  1.77  *10  1 1 ,  and  7.08  •  10  10  F/m, 

res pecti vely 

For  simplicity  it  is  assumed  that  each  of  these  constants  is  uniform 
within  each  region,  and  that  unpaired  charges  occur  solely  at  the  regional 
boundaries . 

Since  this  is  a  dc  analysis,  Faraday's  Law  is  expressed  as: 

V  x  E  =  .  0 

3t 

The  electric  field,  thus,  is  uniquely  represented  by  a  scalar 
potential  that  must  obey  Possoin's  equation: 


V2$  = 


since  there  are  free  charges  only  at  boundaries,  in  each  region  this 
expression  reduces  to  Laplace's  equation: 

V2<jj  -  0 


Membrane  transport  experiments  take  place  over  hours.  Since  the 
charge  relaxation  time  for  these  materials  is  so  fast,  dc  steady  state 
conditions  prevail.  From  charge  conservation: 

V. j  +  -i£L  =  o  =*  V*J 
3t 

Since  charge  transport  is  dominated  by  Ohmic  conduction: 

n*(a1E1  '  a2E2)  “  0 

a  E.  (r  -  RJ  ”  aE  (r  -  R  ) 
e  1  r  e  m  2r  e 

a  E_  (r  -  R  )  -  o  E_  (r  -  R  ) 
m  2r  c  c  3^  c 

and  since  these  fields  are  conservative: 
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♦2(r  -  Rc>  “  'J)3^r  “  Rc) 

These  boundary  conditions  yield  the  following  four  expressions: 
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Therefore: 


Av  = 
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A3  R  (15 oo  ♦  12o2)  +  62R  2(9oo  ♦  l8o2)  ♦  54(6oo  +  3o2)  +  96R,3o2 
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This  can  be  reduced.  A  a  *  Boa  +  Co  *  Ao  if  the  first  term  is  much 

m  m 

Ao2  Ao2 

larger  than  the  other  two,  i.e.  ■= -  >>  1  and  - r  >>  1. 

Boom  r*  2 

m  co_ 

m 

A  °rn  (3  2 

Equivalently,  —  >>  —  and  -7  >>  ( — )  .  Because  of  the  low  membrane 
Bo  A  0 

conductivity  this  is  in  fact  the  case.  Therefore: 
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Likewise,  an  expression  of  the  form  AR  +  BR  6  +  CR  6  +  D6  can  be 

c  c  c 

3  A 

simplified  to  AR  if  its  first  term  is  much  larger  than  the  others,  i.e.  -r 
c  D 

>>  77  >>  (r^-)2,  and  77  >>  (^-)3.  Due  to  the  small  membrane  width  this  is 

R  C  R  DR 
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true  for  cells.  Therefore: 
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Normalizing  each  term  in  this  expression  gives: 
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Thus  if  -z-  >>  —  then 
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Av  =  -6 —  E  cos 9 
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Likewise  the  solution  for  AEc  is: 
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Making  the  same  approximations  gives: 
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Plugging  in  the  numbers  gives: 

AE  «  . 00^8  E  cos 9  for  a  =*10^  mho/m 
com 

-5  -8 

AE  *  4.48  •  10  E  cos 9  for  a  -  1 0  mho/m 
c  o  m 

Therefore,  a  225  mv/m  source  would  induce  an  internal  field  less  than  or 
equal  to  1  mv/m. 

The  expression  for  this  upper  bound  of  electrophoretic  velocity  becomes: 


V  <  6.54  •  1o”12(As2/kg)  E  0030 
e  —  o 


For  a  225  mv/m  source  the  induced  electrophoretic  velocity  would  be: 


Ve  <  1.47  •  10'12  (As2V/kg) 


Ill.g  The  Electrophoretic  Velocity  using  an  AC  Spherical  Cell  Model 


A  slightly  more  accurate  model  can  be  derived  when  the  alternating 
nature  of  the  exciting  fields  is  accounted  for.  As  long  as  the  exciting 
source  remains  within  the  electrostatic  domain,  [15]  if  a  spherical  cell 
were  to  be  exited  by  a  single  frequency  source,  all  of  the  geometric 
boundary  conditions  would  remain  the  same.  Laplace's  solution  would  have 
the  same  form,  but  the  two  time  dependent  conditions  derived  from  charge 
conservation  would  be  modified: 


V.j  ♦  ^L7.D  =  o 


At  the  spherical  boundaries: 
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With  cosinusoidal  input  these  boundary  conditions  become  in  phasor  form: 


(a1  *  joJc1  )Elp  =  (  a 2  +j^2)E2r 
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The  only  difference,  thus,  between  the  ac  solutions  and  the  dc 


counterpart  is  that  a  is  replaced  by  a  +  jwe  .  Their  solution,  therefore, 

X  XX 


is  the  same  as  the  dc  solution  with  a  replaced  with  o  +  jwe  .  The 


approximated  forms  of  these  equations  are: 
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Rearrangement  yields: 
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Since  —  >>  — ,  with  respect  to  frequency,  the  largest  internal  field  and 
0 
m 

the  smallest  membrane  voltage  drop  will  occur  at  the  largest  experimental 
frequency  using  this  model.  Our  experiments  have  ranged  up  to  100Hz.  In 
which  case 

AE  =  .00448  E  cose  for  a  =  1  0  ^  mho/m 
com 

-5  -8 

AE  »  4.48  •  10  E  cos9  for  a  =10  mho/m 
c  o  m 

For  frequencies  at  100Hz  and  below,  both  the  ac  and  dc  expressions  yield 
the  same  numerical  results  for  three  significant  digits.  So  the  same 
electrophoretic  velocity  as  well  would  be  obtained. 

Ill.h  Upperbound  for  Electrophoretic  Velocity  of  a  Cell  of  an  Arbitrary 

Elliptical  Geometry 

When  viewing  transport  data,  although  instructive,  a  spherical  model 
for  a  cell  is  clearly  not  accurate.  Fibroblasts  growing  in  monolayer  lie 
very  close  to  the  glass  plane,  while  they  can  reach  up  to  2.5mm  in  the 
other  directions.  From  observations  it  is  clear  that  an  ellipsoidal  model 
would  be  more  appropriate.  But  these  generalized  spheroids  do  not  allow  a 
separation  of  variables  for  Laplace's  equation,  and  therefore  they  can  not 
be  solved  analytically.  [8]  On  the  other  hand  they  can  not  simply  be 
ignored  since  most  data  indicates  that  for  a  given  length  in  the  direction 
of  the  electric  field  the  intracellular  field  is  lowest  with  a  sphere.  [16] 
So  the  electrophoretic  effects  that  are  not  possible  in  the  sphere  could 
conceivably  occur  with  a  more  oblate  or  prolate  model.  We  have,  thus, 
developed  an  expression  for  the  upper  bound  for  an  intracellular  electric 
field  in  a  convex  geometry.  Note:  this  expression  should  only  be  used  as  an 
upper  bound.  Due  to  the  difficulties  with  the  geometry,  it  is  undoubtedly 
considerably  higher  than  the  real  field.  But  such  a  calculation  has  its 
purposes . 

Intracellular  electric  fields  for  generalized  geometries  can  not  be 
solved  exactly,  but  if  the  voltage  drop  at  all  points  on  the  cell  membrane 
is  known,  the  current  through  any  region  of  the  membrane  would  be  known: 


v  (u  ,u  )a 
m  1  2  m 


du1  du^ 


A  closed  contour  at  a  constant  value  for  one  of  the  ellipsoid's  principal 
axes,  call  it  A,  is  a  two  dimensional  shell  (ellipse  or  a  circle  depending 
on  the  axis).  By  charge  conservation,  the  total  amount  of  current  flowing 
through  the  membrane  on  either  side  of  A  must  equal  the  total  amount  of 
current  flowing  the  planar  surface  bounded  by  A. 


where  i  =  the  current  flowing  through  the  planar  surface  bounded  by  A,  ift, 
=  the  current  flowing  through  the  smaller  region  of  the  membrane  bounded  by 


bounded  by  A.  Using  Ohm's  Law  and  uniform  conductivity: 
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where  E&vA  is  the  average  value  of  the  normal  component  to  the  planar 
surface  of  A  of  the  intracellular  electric  field,  and  S_  is  the  Surface 

cL 

area  of  the  planar  region  bounded  by  A. 
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Where  S  is  the  surface  area  of  the  smaller  portion  of  membrane  bounded  by 
the  contour  A. 
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So  for  any  given  planar  contour  there  is  a  bound,  D,  for  the  maximum 

component  of  the  electric  field  perpendi cular  to  the  plane.  To  find  the 

bound  for  the  maximum  component  of  an  electric  field  within  the  ceil,  the 

contour  in  which  D  is  maximized  must  be  found.  That  contour  is  the  contour 

which  maximizes  S  /S  . 

m  a 


III.h.1  An  Upper  3ound  for  3  /S 

m  a 

As  shown  in  Figure  9,  a  general  spheroidal  geometry  can  be  expressed 
as : 


1 


3y  symmetry  the  bound  f or  a  contour  at  x=f  will  equal  the  bound  for  a 
contour  at  y»f.  So  there  are  two  indeoendent  regions  where  this  problem 
must  be  solved  (Figure  10).  First  for  z=c  c>0  -  z§-c:  the  above  expression 
yields  a  contour  a  which  is  a  circle  with  radius: 
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Call  this  bound,  F.  F's  extremum  with  respect  to  c  can  be  found  by  taking 
the  first  derivative. 

3F  -2b2  +2bc 

3c  =  2  c2  1 

bad-  — p)2 

b 

Thus,  F's  extremum  comes  at  c  =  b.  The  second  derivative  is  not  finite  at 
this  point,  but  since  this  is  the  only  extremum  of  F,  F  is  monotonic  on 
both  sides  of  c=b .  The  two  constrained  extrema  for  the  interval  0  _<  c  £  b 
will  come  at  c  =  0  and  c  -  b .  One  of  which  necessarily  will  be  the 
constrained  maximum  and  the  other  will  be  the  constrained  minimum. 


F(c  =0)  =1  +  — 
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F(  c  =  b )  =  1  + 


_ 2(b-c) 

2  2  — 
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Applying  L'  Hopital's  Rule  gives: 


F(c  -  b)  -  1 
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III.h.2  An  Upperbound  for  the  Surface  Area  Ratio  in  the  Other  Case 


If  x  -  f: 
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The  area  of  the  membrane  is  less  than  that  of  the  ellipse  +  the 
circumference  times  the  height  a-f. 

S  <  S  +  C  (a-f) 
m  a  a 

The  circumference  of  the  ellipse  is  less  than  the  circumference  of  a 

rectangle  of  sides  2y__  and  2z_„. 
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Call  this  bound  G.  Taking  its  first  derivative  gives: 


9G  Ma(f-a) (a+b) 

3f  "  o  ?  o  3 

na^bd  -  f  /a“)2 

The  extremum  of  G  comes  at  f  =  a.  Unfortunately  the  second  derivative 
is  discontinuous  at  this  point.  But  f=a  is  the  only  extremum  of  G; 
therefore,  G  is  monotonic  on  both  sides  of  f-a.  In  the  range  of  interest  0 
£  f  £  a,  the  constrained  extrema  will  thus  occur  at  f=0  and  f=a.  One  of 
these  will  be  the  constrained  maxima;  the  other  the  constrained  minimum. 
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Applying  L'Hopital's  Rule  gives: 


G(f-a)  =  1 
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III.h.3  The  Cellular  Configuration  that  Yields  the  Largest 

Induced  Membrane  Potential 

A  separation  of  variables  is  not  possible  for  general  elliptical 
geometries,  so  an  exact  analytic  solution  is  not  currently  possible.  A 
technique  for  approximating  the  potential  that  is  sometimes  used  is  to 
assume  that  the  membrane's  conductivity  is  zero.  This  provides  boundary 
conditions  which  are  separable  in  spheroidal  coordinates.  [8] 

There  is  justification  for  assuming  negligible  conductivity.  In 

_6  _  Q 

general  the  membrane  conductivity  might  be  10  -  10  mho/m,  whereas  the 

cytoplasmic  and  media's  conductivity  are  at  Imho/m. 

For  this  particular  application  there  might  even  be  more  of  a  reason 

for  making  such  a  calculation  for  V  .  When  Bernhardt  and  Pauly  used  the 

max  J 

above  approximation  they  found  that  the  maximum  voltage  drops  occur  when 
the  cells  have  a  line  of  symmetry  parallel  to  the  electric  field.  And  they 
found  that  this  maximum  occurs  at  the  line  of  symmetry.  Let's  look  at  an 
electric  field  line  at  this  line  of  symmetry  with  an  area  of  dA.  From 
steady  state  solutions  it  is  known  that: 

(c  +  jwe  )E  *dA  -  ( o  +  jwe  )E  *dA  ■  (a  *juie  )E  *dA 
e  ee  m  mm  c  cc 

Therefore  as  a  is  decreased,  E  would  be  expected  to  increase.  It 
m  m  r 


would  not  be  expected  to  Increase  by  the  same  proportion  as  the  decrease, 
since  the  other  parameters  would  also  be  affected.  But  an  increase  would 


definitely  be  expected.  Since  V  »  E  6  ,  V  would  also  be  expected  to 

ra  m  m  m 


increase.  Therefore,  as  a  -*■  0  we  would  expect  that  V  might  reach  a 

m  m 


maximum . 


Intuition  aside,  numerical  data  for  several  geometries  indicates  that 
these  approximations  are  reasonable  for  intracellular  and  extracellular 
electric  fields.  [16]  Thus  one  would  expect  it  to  be  valid  for  the  voltage 
across  the  membrane  as  well. 


Since  it  appears  the  maximum  voltage  drop  occurs  when  an  axis  of 
symmetry  is  parallel  to  the  imposed  electric  field  [8],  the  data  from 
Bernhardt  and  Pauly  can  be  used.  For  a  given  height  in  the  direction  of  the 
electric  field,  the  highest  membrane  voltage  drop  occurred  with  an  oblate 
spheroid  with  its  small  axis  in  the  direction  of  the  electric  field. 


For  cells  in  monolayer,  however,  this  is  a  ridiculous  model.  The 
shortest  axis  is  perpendicular  to  the  cover  glass  which  is  also 
perpendicular  to  the  electric  field.  The  largest  drop  for  a  cell  in  a 
reasonable  geometry,  therefore,  is  that  of  a  cylinder  with  its  large  axis 
perpendicular  to  the  electric  field.  Bernhardt  and  Pauly  found  this  to  be 
the  second  largest  drop. 


III.h.4  The  Maximum  Change  in  Membrane  Potential 


Using  Laplace's  equation  and  making  the  same  assumptions  as  for  the  ac 
spherical  case  we  get  the  following  solutions  to  the  electric  field  being 
applied  to  a  cell  with  a  cylindrical  geometry: 
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With  boundary  conditions  of: 
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Thus  there  are  four  linear  equations  for  the  four  unknown  parameters: 
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Solving  these  equations  simultaneously  and  plugging  in  numbers  gives: 

AV  -  .002m  •  E 
m  o 

III.h.5  The  Upper  Bound  for  the  Intracellular  Electric  Field 

Accumulating  the  results  of  the  last  few  sections,  the  upper  bound  for 
the  average  value  of  any  coordinate  component  of  the  electric  field  is: 


am  Vmax  Sm  _ ,.  .  2b  ,  .  2(a  +  b), 
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IV.  Distortion  and  Rectification  Due  to  Electrode-Electrolyte 

Int er act ions 


IV. a  Introduction 

Relating  the  applied  electric  field  that  the  cells  perceive  to  the 
current  and  voltage  that  is  applied  at  the  electrodes  is  very  important  to 
interpreting  the  results  of  our  proline  incorporation  data.  For  this  reason 
we  have  undertaken  a  study  that  analyzes  the  amount  of  distortion  and 
rectification  that  occurs  in  an  applied  field  with  respect  to  its  source. 

IV. b  Background 

In  the  proline  experiments  an  electric  field  was  applied  across  a 
region  of  a  chamber  by  two  electrodes  placed  on  either  side  of  the  region 
(Figure  11).  For  this  regional  field  to  be  a  pure  sinusoid,  devoid  of 
harmonic  distortion  or  rectification,  the  applied  voltage  and  current  must 
belong  to  particular  families  of  solutions.  Namely,  the  applied  voltage 
must  possess  a  certain  normalized  spectrum,  and  the  applied  current  must 
possess  a  particular  normalized  spectrum  that  are  solutions  to  the 
constraints  imposed  by  the  chamber. 

•  Current  in  a  region  is  related  to  the  voltage  across  the  region  by 
the  expression: 


Figure  11 


Electrode  area  =  1.22cm-2. 


Figure  12 


I  -  Y(V)V 


where  Y  is  the  admittance  of  the  region.  By  using  admittances,  charge 
conservation  can  be  expressed  as: 


1 1  -  0 

for  currents  entering  a  node.  From  Figure  11  then  it  can  be  seen  that: 

^applied  *bath 

the  media  in  the  region  can  be  accurately  modeled  as  having  a  real, 
constant  conductivity.  Then  assuming  the  region's  geometry  results  in  a 
uniform  field: 
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The  ratio  of  the  region's  field  to  the  applied  current  is  constant; 
therefore,  they  possess  the  same  normalized  spectrum.  For  a  purely 
sinusoidal  field,  the  applied  current  must  be  purely  sinusoidal. 

IV. c  Voltage 


Ampere's  and  Faraday's  Laws  for  this  electro-quasistatic  system  can 
be  expressed  as: 
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for  a  uniform  field: 
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Expressions  such  as  the  Butler  -  Volmer  Equation  and  experimental 

data  indicate  that  Y  (E  )  is  not  a  constant.  Therefore,  the  ratio  is  not 

e  e 

constant  with  respect  to  E.  The  region's  electric  field  and  the  applied 
voltage  do  not  share  the  same  normalized  spectrum.  A  purely  sinusoidal 
applied  voltage  will  not  produce  a  purely  sinusoidal  E;  otherwise,  they 


would  share  the  same  normalized  spectrun.  sk 

To  get  a  purely  sinusoidal  regional  field,  one  must  apply  a  purely 
sinusoidal  current  and  a  more  complicated  voltage  —  a  voltage  that  has 
frequency  components  other  than  the  fundamental.  It  can  be  said  then  that 
the  applied  voltage  is  harmonically  distorted.  Since  it  does  not  possess 
the  same  spectrum  as  a  pure  sinusoid,  it  is  possible  that  the  applied 
voltage  might  also  be  rectified  with  respect  to  the  regional  field. 

IV. d  Rectification 

If  the  regional  electric  field  is  not  rectified  with  respect  to  the 
applied  voltage,  then  in  all  but  the  most  pathological  cases,  the  regional 
field  must  be  an  odd  function  of  the  applied  voltage;  that  is 

Ec(‘  Vapplied)  *  "  Ec( Vapplied) ' 

Since 

I  =  E  o  A*-  I 
c  c  c 


this  condition  can  equivalently  be  expressed  that  I  must  be  an  odd 


'unction  of  V 
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Since  the  admittance  is  a  one-to-one  mapping  of  voltage  to  current, 
Y(V)  can  be  expressed  as  Y(I)  and  Z(I)  can  be  calculated. 

V  „nH^U)  -  Z  (1)1  ♦  Z-(I)I  ♦  Z  I 

applied  el  e2  c 

V  i.JD  -  (Z.(I)  ♦  Z  _(I)  ♦  Z  )I 

applied  el  e2  c 

The  sum  Z  .  (I)  +  Z  _(I)  needs  to  be  an  even  function  of  I  to  avoid 

ei  e2 

rectification. 

From  Figure  11  the  current  I’s  positive  direction  is  defined  to  be  the 
path  within  the  region  from  electrode  1  to  electrode  2.  Each  electrode  can 
be  described  as  having  a  characteristic  impedance,  zChar^ Ichar ^ ’  for 

current  that  is  defined  as  positive  when  flowing  from  the  electrode  to  the 
electrolyte.  For  electrode  1,  IQhar  ■  I.  For  electrode  2,  Ichar  *  -I. 


Therefore; 


ze,(D  .  ze2(i)  -  Z0har.8,c»  *  zonar.e2(-i) 


Thus,  if  the  characteristic  impedance  functions  of  the  two  electrodes 
are  the  same,  this  sum  will  be  an  even  function  of  I  —  no  rectification 
will  occur.  If  they  are  not  the  same,  then  in  all  but  extremely 


pathological  cases,  there  will  be  rectification. 


Real  electrodes  are  not  likely  to  have  identical  impedances. 
Rectification  of  the  regional  field  with  respect  to  the  applied  voltage  is 
expected. 

To  produce  a  purely  sinusoidal  electric  field  in  the  chamber's 
region,  a  purely  sinusoidal  current  and  voltage  that  is  rectified  and  is 
not  purely  sinusoidal  must  be  applied.  Under  these  circumstances,  a  current 
source  is  much  better  to  use  than  a  voltage  source.  A  purely  sinusoidal 
current  gives  a  purely  sinusoidal  field;  a  purely  sinusoidal  voltage  gives 
a  harmonically  distorted  and  rectified  field. 

IV. e  Experiments 

Due  to  the  previous  analysis,  distortion  and  rectification  for  an 
applied  sinusoidal  voltage  is  expected.  As  a  result  of  the  imperfect  nature 
of  current  sources,  the  distortion  and  rectification  of  its  normalized 
spectrum  is  expected.  However,  these  quantities  are  expected  to  be  less  for 
the  current  source.  We  have  conducted  experiments  to  quantify  the 
distortion  and  rectification  of  the  regional  electric  field  for  real, 
sinusoidal  current  and  voltage  sources. 

The  chamber  used  is  shown  in  Figure  12.  It  consists  of  a  cut-open  10  cc 
syringe  that  is  closed  by  two  rubber  stoppers  at  the  opposing  ends.  Each 
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rubber  stopper  has  a  circular  1.22  cm  electrode  attached  to  it,  and  each 
stopper  was  movable  along  the  main  axis  of  the  syringe. 

The  electrodes  were  made  of  bright  platinum  metal  that  was  buffered 
with  an  emory  cloth.  They  were  affixed  to  the  stoppers  with  Epoxy.  After 
which,  they  were  repeatedly  washed  with  a  Q-tip  soaked  in  IN  HCl/HNO^  and 

wiped  with  an  emory  cloth.  The  emory  cloth  increased  the  effective  surface 
area  of  the  electrode  and  removed  most  of  the  epoxy  that  might  have 
attached  to  the  electrodes. 


Media  was  introduced  into  the  syringe  through  a  small  hole  that  was 
made  in  the  side  of  the  syringe.  Each  day  media  was  loaded  and  removed 
through  the  hole. 


The 

experimental 

configuration 

is  shown 

in 

F i gur  e  13.  The 

test 

resistors 

were  chosen 

so  that  the 

magni tude 

of 

the 

impedance  of 

the 

resistor 

corres ponded 

to  ten  to  fifty  percent 

of 

the 

magnitude  of 

the 

path's  total  impedance.  From  the  previous  analysis,  it  is  known  that  the 
current  through  the  resistor  is  linearly  related  to  the  electric  field  in 
the  media.  Since  current  is  linearly  related  to  voltage  through  a  resistor, 
it  is  known  that  the  voltage  across  the  resistor  is  linearly  related 
without  rectification  to  the  electric  field  in  the  media.  Thi3  voltage  then 
possesses  the  same  normalized  spectrum  as  the  electric  field  in  the  media 
—  from  it  the  harmonic  distortion  and  rectification  of  the  media's 
electric  field  with  respect  to  a  pure  sinusoid  was  determined. 

The  voltage  drop  was  converted  to  a  digital  signal  and  stored  by  using 
an  Analog  to  Digital  Converter  interfaced  to  an  IBM- XT  computer.  The 
frequency  spectrum  was  calculated  using  the  Asyst  scientific  software 
package . 


IV. f  Data 


Experimental  runs  were  conducted  using  two  different  media:  . 15M  KC1 
and  Dulbecco's  Modified  Eagle  Media  (Gibco  —  IX  preparation).  The  voltage 
source  was  a  Wavetek  waveform  generator,  and  the  current  source  was  the 
custom-built  current  source  that  was  used  in  our  proline  incorporation 
experiments.  Both  sources  are  digital  generators  with  high  resolution  and 
stable  frequency  outputs.  The  Wavetek  generator's  output  was  found  to  not 
be  purely  sinusoidal.  When  the  chamber  was  replaced  by  a  resistor  in  the 
configuration  of  Figure  13.  a  0.2%  dc  offset  and  a  harmonic  distortion  of 
less  than  2%  were  found  at  each  of  the  frequencies  and  current  densities 
for  which  the  experiment  was  run.  This  "background"  distortion  and 
rectification  was  subtracted  from  each  of  the  stored  spectra  to  yield 
"corrected"  spectra. 

The  following  frequency  sinusoids  were  applied:  .1,  1,  10,  and  100 
Hz.  The  following  peak  current  densities  were  applied  at  each  frequency  for 
both  sources:  .001,  .01,  .1,  1,  10  mA/cm2.  From  the  resulting  spectra  the 
percent  harmonic  distortion  and  rectification  were  calculated.  Percent 
harmonic  distortion  was  defined  as  one  hundred  times  the  magnitude  of  the 
dc  component  of  the  corrected  and  normalized  spectrum  (normalized  with 
respect  to  magnitude  not  intensity).  Percentage  harmonic  distortion  was 
defined  as  one  hundred  times  the  sum  of  the  magnitude  of  the 
non-fundamental  harmonics  of  the  corrected  and  normalized  spectrum. 

The  results  are  plotted  in  Figures  14-15.  Percent  rectification  was 
found  to  decrease  with  increasing  current  densities  and  with  increasing 
frequencies.  Percent  distortion  increased  with  increasing  current 
densities,  but  decreased  with  increasing  frequency. 

This  raised  a  major  question:  what  is  the  mechanism  of  the  frequency 
dependence  of  the  metal-electrode-electrolyte  rectification?  To  explain 
this,  future  work  will  quantitatively  examine  the  specific  electrochemical 
reactions  that  have  rate  constants  on  the  order  of  .1  to  10  seconds. 
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